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1. INTRODUCTION 

Protection engineers face challenges in certain applications such as islanding detection in a 
distribution network. Island is a situation where a part of a utility grid is energized by DGs and electrically 
separated from rest of the network. It is very much important to detect islanding condition correctly & 
quickly, otherwise risk of damage to the power system rises & safety hazard for the personnel becomes a 
matter of concern. Several power systems face the issue of low frequency oscillation due heavy load 
condition or to system interconnection. The stability of a power system is a non-linear, dynamic 
phenomenon, and it depends on poorly damped low frequency oscillations. These oscillations play vital role 
in the analysis of stability of the system. If these oscillations are not damped sufficiently, an unstable 
operation may occur and it may lead to a network collapse. Therefore, it is essential to monitor the modal 
parameters of the low frequency oscillatory signals for dynamic system security [1]. 

Reference [2] presents several islanding detection techniques such as active inverter-resident 
techniques, passive inverter-resident techniques, communications-based techniques, and utility level 
techniques for distributed power generation systems. In Reference [3], an islanding scheme based on wavelet 
singular entropy is implemented in a microgrid with DG. A hybrid analyzing method which analyzes the d- 
axis voltage component in 2 ways for detecting the islanding precisely is presented in [4]. A hybrid islanding 
detection method for microgrids (MGs) with various connection points to smart grids (SGs) is proposed in 
[5], and it based on probability of islanding calculated at the SG side and sent to the central control for MG. 
Reference [6] analyzes the sensitivity of 16 power system parameters which are used in passive methods to 
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detect non-islanding and islanding events. An islanding detection method based on islanding discrimination 
factor is proposed in [7], and it is derived from the superimposed components of voltages. An islanding 
detection method based on average absolute frequency deviation value is presented in Reference [8]. 

An islanding detection method for MGs considering small scale synchronous generators is presented 
in [9]. Tuned filters connected at the DG terminals are utilized for islanding detection in MG is proposed in 
[10]. An islanding detection method based on modal components of voltages is discussed in [11]. An active 
islanding detection method for an inverter-based DG is proposed in [12]. A passive local multi-criteria based 
islanding detection technique with fault detection as well as islanding verification logic is proposed in [13]. 
Reference [14] presents an efficient approach for building decision trees-based intelligent relay. A new large 
performance islanding search sequence method implemented to 4 islanding detection approaches is proposed 
in [15]. A new passive-based anti-islanding technique for both synchronous machine & inverter based DGs is 
proposed in [16]. 

It can be noticed from the literature survey that there is a pressing requirement for fast and accurate 
algorithm for islanding detection. Therefore, in this paper, two signal processing based islanding detection 
methods, i.e., DWTwith ANN, and S-transform with ANN are proposed. Here, PSO based feature selection 
(FS) algorithm is used. The suitability and effectiveness of proposed islanding detection methods has been 
examined on distribution network with DGs. 


2. PROPOSED SIGNAL PROCESSING TECHNIQUES FOR ISLANDING DETECTION 
In this paper, DWT and S-transforms are applied for detecting the islanding in the 
distribution network. 


2.1. Islanding detection using discrete wavelet transform (DWT) with ANN 

Wavelet transform (WT) is an adaptive signal processing technique for non-stationary signals and it 
extractsinformation from the data series. Choice of desired mother wavelet is vital for good performance, and 
Daubechies is considered as the appropriate one for analyzing the transient event [17]. Here, WT is applied 
for islanding detection. Signal is collected from the relaying end of distribution line is decomposed into 8 
level by WT, and it is depicted in Figure 1. Sampling frequency is considered as 30 kHz, therefore, 8 level 
decomposition can be performed [18]. For the accuracy of the model, sampling frequency of 30 kHz is 
required. Decomposition process depends on sampling frequency (Fs) and number of samples 

Decomposition process with digital filtering approach is depicted in Fig.1. In this figure, X[n] is 
signal, g[n] is high pass filter, and h[n] is low pass filter [19-20]. ‘dl’and ‘al’ are 1* level detail and 
approximation coefficients, respectively. 


Level-1 Level-2 Level-8 d8 


0.117 ~ 0.234 
kHz a8 





0 ~ 15 kHz 0 -7.5 kHz 0 ~ 0.117 kHz 


Figure 1. Digital filtering technique of wavelet transform (WT). 


ANN is a computational model that simulates the structural and functional aspects of biological 
neural network [21]. In this work, a multilayer feed-forward neural network (FFNN) with back propagation 
training algorithm is applied. Flow chart of islanding detection technique in a distribution network using WT 
and ANN is depicted in Figure 2. At the targeted distribution generation (DG) location, the negative sequence 
current/voltage signals are retrieved. These signals are decomposed using DWT [22]. The obtained statistical 
features from the signals are reconstructed from detailed coefficients. Feature set with best predictive power 
is chosen by using PSO algorithm. Thereafter, the train & test data is generated from wide variation of 
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loading condition and it is normalized [23]. ANN is trained with the selected features as input. Then, the 
testing is performed with the trained neural network. 
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Figure 2. Islanding detection technique using DWT and ANN. 


2.2. Islanding detection using S-Transform with ANN 

S-transform represents the time-frequency relationship with instantaneous values of frequency, 
phase and amplitude [24]. It gives the absolutely referenced phase information & frequency invariant 
amplitude response & provides better signal clarity for transient signal. Generates a contour which is simple 
to visualize, whereas WT requires standard multi-resolution analysis. However, it requires higher complex 
computation [25]. S-transform confines the real & imaginary components, phase as well as amplitude 
spectrums independently. It is termed as absolutely referenced phase information. S-transform is used for 
several applications such as fault classification, location and modal analysis of signal. 

Figure 3 depicts the flow chart of proposed islanding detection scheme using S-transform & ANN. 
Negative sequence voltage/current signals are acquired at the targeted DG location. S-transform is applied to 
these signals, and energy is determined [26-27]. The cumulative sum (CUMSUM) of energy signal is 
checked for the convergence condition, and this CUMSUM is used for islanding detection. 
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Figure 3. Islanding detection technique using S-transform and ANN. 


2.3. Feature extraction 

It is used for reducing data matrix size by transforming into features. From decomposed coefficients 
of DWT, six statistical features are extracted. For WT, the number of reconstructed coefficients are 8 (i.e., 
d1-d8) and the total feature matrix is 96 (i.e., 8'decomposition level x 2nos of signals x 6 nos of features). 
Total feature set for each signal for S-transforms is 18 features [28]. In S-transform, 3 parameters, i.e., 
frequency, phase and amplitude are considered. 

The transient energy signal has large value as compared to normal signal. Standard deviation of an 
undistorted signal is equal to one, however, for transient signal, this value deviates from one. During the 
islanding event, the signal mean becomes non-zero, whereas its value is zero for undistorted signal [29-30]. 
For an undistorted signal, the value of kurtosis is 3, whereas its value is greater than 3 for a transient signal. 
For an undistorted signal, the value of skewness is zero & non-zero for distorted case. The value of entropy is 
large value for transient signal, and it is less for an undistorted signal. 
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3. FEATURE SELECTION 

Feature selection (FS) is selecting small sets of features which determine’s target properly. It 
reduces the complexity of prediction & learning technique & makes prediction accuracy to rise. In various 
supervised learning activities, most of the features shows redundancy [31]. Therefore, FS is a significant 
approach towards learning of large multi-featured data matrix. 


3.1. Feature selection using PSO algorithm 

PSO is stochastic optimization approach motivated by the behavior of bird flocking or fish 
schooling. In PSO, the potential solution is called particles, which interacts among themselves to find global 
optimal solution. In every iteration, each particle has its own position & velocity, and is updated by two best 
values called pbest&gbest. For further details of PSO algorithm, the reader may refer reference [32]. The 
flow chart of PSO based FS is depicted in Fig. 4. 
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Figure 4: Flow chart of PSO based feature selection. 





3.2. Parameter setting for generating train & test data matrix 

The train & test data matrix is normalized between [0,1]. For each signal, 48 features (8 level 
decomposition x 6 features) for DWT, and 18 features (3 modal parameters X 6 features) for S-transforms 
are considered. Optimal features are obtained from feature selection method, and they are 8 for DWT and 4 
for S-transform. In this paper, two types of training data sets are considered. The first data set is formed by 
the combination of 7 resistances (0Q, 19, 5Q, 10Q, 25Q, 35Q, 45Q) and 6 inception angles (10°, 20°, 30°, 
40°, 50°, 85°). Here, ten forms of event are considered with 300 locations [33]. Therefore, the total of 126000 
(i.e., 10 X 300 X 42) training data matrix is formed. The second data set is formed by the combination of 9 
resistances (2Q, 4Q, 6Q, 9Q, 12, 20Q, 30, 40Q, 50Q) and 8 inception angles (5°, 11°, 17°, 24°, 32’, 45°, 
60°, 90°). Therefore, the total of 216000 (i.e., 10 x 300 x 72) training data set is generated. In this paper, the 
accuracy (in %) can be calculated by using [34], 


|Total number of samples — Misclassified| 


0 ee 
ENE Total number of samples ree a) 


4. SYSTEM DESCRIPTION 

The system under study for islanding detection in a distribution network is depicted in Figure 5. The 
grid data considered in this work has a frequency of 50 Hz, voltage of 120 kV, zero sequence parameters: Ro 
is 1.728 Q/km and Lo is 0.055 H/km, and the positive sequence parameters: R: is 0.576 Q/km and L; is 0.018 
H/km. Distribution line data considered in this work has the line length of 30 km, zero sequence parameters: 
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Ro is 0.826 Q/km, Lo is 6.64 mH/km, Co is 0.1002 nF/km, positive sequence parameters: Rj is 0.2306 Q/km, 
Lı is 2.1 mH/km and C; is 0.2266 nF/km. The DG parameters has generator data: wind speedis10 m/s, 
number of wind turbines are 6, nominal power is 9 MW (i.e., 6 x 1.5 MW), voltage is 400 V, frequency is 50 
Hz; DC bus voltage regulator gain has kọ is 8, k; is 400; grid-side converter current regulator gain has k, is 
0.83, k; is 5; speed regulator gain has k,, is 3, k; is 0.6; rotor side converter current regulator gain has ky is 
0.6, k; is 8; reactive power & voltage regulator gain are 0.05, 20; pitch controller gain kpis 150; converter 
data has nominal DC bus voltage has 1150 V and DC bus capacitor is 10000 uF. 


9 MW DFIG BASED DG-1 
WIND FARM 


L2 
25 kV/ 400 V TR-2 





TR-1 






L1 25 kV/ 400 V TR-3 


L3 


9 MW DFIG BASED 
WIND FARM DGZ 


Figure 5: System under study for islanding detection. 


From figure 5, it can be observed that the system under study has two 9 MW wind farms driven by 
wind turbine. Each 9 MW wind farm consists of six wind turbines of 1.5MW capacity connected to 120 kV 
grid through 25kV, 30 km feeder. The sampling frequency considered is 200 kHz &system frequency is 50 
Hz. Therefore, there is 4000 samples/cycle. Load demands are varied at DG end as well as at point of 
common coupling. Current samples are retrieved at DG-1 and DG-2 ends. In this work, the wind speed is 
considered as 10 m/s. Here, two cycles of current signal is considered, one just before islanding and another 
after islanding. 


5. RESULTS AND DISCUSSION 
In this paper, 18 different combinations of loading conditions are utilized to generate train & test 

data. For each signal, 48 features (8 level decomposition x 6 features) are considered for DWT. Ten optimal 
features are selected from the PSO based optimal feature selection method, and they are fed to ANN for 
training. This trained neural network is used for testing purpose. To avoid over fitting of data, ten fold cross 
validation is carried out [35]-[36].70% training data and 30% testing data will be selected randomly for the 
test set. Because of this, some observations may not be selected in the test set, whereas others may be 
selected more than once. This results in test set overlapping [37]. To overcome this situation, cross-validation 
is performed. In this paper, the cross validation is performed, and checked that for all the test cases error is 
within the limited range. Also, the proposed method in the paper is more robust than cross-validation, 
because in cross-validation the operating condition is already been seen by the neural network, whereas in the 
proposed technique, operating condition for test data is altered from trained one [38]. As mentioned earlier, 
in this work, two islanding detection approaches are used, and they are 

e Method 1: Islanding detection using DWT with ANN. 

e Method 2: Islanding detection using S-Transform with ANN. 
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5.1. Islanding detection by using method 1 
As mentioned earlier, in this method, DWT with ANN is used for islanding detection in distribution 
network. Table 1 presents proposed methods accuracy for islanding & non-islanding conditions. 


Table 1. Accuracy of proposed approaches for islanding & non-islanding conditions 








Condition No. of samples Correct identification Mis-identification Accuracy (%) 
Method 1 Method2 Method 1 Method 2 Method 1 Method2 Method 1 Method 2 
Islanding 18 18 18 17 0 1 100 94.4 
Non-islanding 54 54 54 49 0 5 100 91 





The obtained output after using discrete wavelet transform (DWT) is fed to an ANN, and it is 
depicted in Figure 6. 


+ Actual 
O ANN output 


41- Island, 0- nonisland 





Number of samples ` 


Figure 6. Islanding detection by using DWT and ANN. 


Figure 6 shows '0' for fault free condition in the system, i.e., under normal operating conditions, and 
it shows '1' when the fault occurs under islanding condition. 


5.2. Islanding detection by using method 2 
Figure 7 depicts the plot for retrieved signal samples taken at the standard frequency of 50 Hz, and 
the fault is occurred after 4000 samples. 


Frequency(Hz) 





1000 2000 3000 4000 5000 6000 7000 8000 
Number of samples 


Figure 7. Islanding detection using S-transform. 


From figure 7, it can be observed that after 4000 samples, there is the detection of islanding 
condition, and hence there is a sudden surge in the system frequency. 
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6. CONCLUSIONS 

This paper proposes an approach for detecting islanding in distribution system. Two islanding 
detection techniques are proposed in this paper based on DWTwith ANN, and S-transform with ANN. From 
the simulation results on considered system for islanding detection in a distribution network shows that the 
DWT in combination with ANN gave 100% accuracy, and it is robust and more accurate than the other 
techniques presented in the literature. In this paper, feature selection is proposed by using PSO algorithm. 
Feature selection makes proposed approach more superior than other methods reported in the literature. 
Extending the present work to meet the problem caused due to the sudden change in load which may create 
false alarm in the islanding detection process is a future research scope. 
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